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d 要 : 在 各 向 腊 性 网 格 下 ， 讨 论 了 一 类 非 线 性 双 曲 型 积分 微分 方程 的 一 个 和 珑 形 非 协调 有 限 元 方法 逼近 ， 
给 出 了 了 半 痪 散 格式 下 的 有 限 元 解 的 收敛 性 分 析 和 误差 估计 。 在 精确 解 适当 光滑 的 前 提 下 ， 利 用 新 
的 技 与 和 精细 估计 得 到 了 其 超 通 近 性 质 。 同 时 利用 插值 后 处 理 技术 导出 了 整体 超收 敛 结果 。 本 文 
的 结论 表明 传统 有 限 元 分 析 中 对 网 格 的 正则 性 要 求 和 对 Ritz-Volterra 投 影 的 依赖 不 是 必要 的 ， 
从 而 进一步 扩展 了 非 协调 有 限 元 方法 的 应 用 范围 。 
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1 引言 


非 线 性 双 曲 型 方程 及 非 线 性 双 曲 想 积 分 微分 方程 这 类 问题 的 研究 在 核反应 动力 学 、 粘 弹性 
力学 、 生 物力 学 、 松 散 介质 中 的 压力 等 实际 问题 的 理论 和 应 用 方面 均 有 一 定 的 价值 。 关 于 这 些 
非 线性 方程 的 有 限 元 方法 的 研究 已 有 许多 上 3 ， 但 大 部 分 要 求 剖 分 满足 正则 性 (B) E < c， 这 
里 hg 为 剖 分 单元 KK 的 直径 ，pk 为 天 的 最 大 内 切 圆 直径 ，c 是 一 个 常数 )。 但 最 近 的 一 些 研究 
成 果 回 表明 ， 这 种 假设 对 一 些 单元 来 说 是 不 必要 的 ， 对 那些 定义 在 窗 边 区 域 上 的 问题 如 果 用 正 
则 性 剖 分 ， 计 算 量 将 非常 大 而 无 法 承受 ， 这 时 若 利 用 各 向 异性 网 格 则 可 以 节约 大 量 的 自由 度 而 
获得 同样 的 收敛 效果 。 另 一 方面 以 往 对 非 线 性 双 曲 型 方程 、 非 线性 双 曲 型 积分 微分 方程 及 一 般 
的 发 展 型 方程 的 研究 几乎 都 是 关于 协调 有 限 元 ， 而 且 Ritz-Volterra 投影 等 传统 方法 是 必 不 可 少 
和 的， 其 误差 估计 也 主要 集中 于 收敛 性 分 析 。 对 于 非 线 性 双 曲 型 积分 微分 方程 如 何 进 一 步 提高 有 
限 元 解 的 精度 ， 仍 是 有 待 解决 的 问题 。 

本 文 把 文献 外 中 的 非 协调 元 应 用 到 非 线 性 双 曲 型 积分 微分 方程 ， 利 用 有 限 元 插值 性 质 讨论 
了 半 离 散 格式 下 非 协 调 有 限 元 解 的 收敛 性 。 同 时 利用 一 些 新 的 技巧 和 精细 估计 得 到 了 超 示 近 性 
质 。 最 后 ， 通 过 构造 适当 的 各 向 异性 播 值 后 处 理 算 子 ， 导 出 了 整体 超收 敛 的 结果 。 


2 单元 构造 
(835 à K = [一 1,1] x [1,1] 上 构造 有 限 元 (K, P, È) 如 下 
X = (61, 65, 65, 04, 05], P = span{1, £, 7] v(£), e(n). 
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其 中 1 1 1 
ô; = — | ôd, i=1,2,3,4, = — | ddédn, t) = =(3t2 — 1). 
6 i; / 9 i 3 ds a; tm e(t) = 3( ) 
容易 验证 ， 对 任意 的 六 € H1( 请 )， 其 插值 函数 可 表示 为 
1 1 、 
Iô = v + T — V4)€ 十 T 一 $1) 十 3 ( T 44-— 2ós)v(£) T 3 (5 +i- 20s)v(1), 
且 具 有 各 向 异性 特征 回 。 
有 限 元 空间 V^ xg SC 
vV’? = Io =wvlkoFk e P, VK €T, | [lv*]Jds 20, F C ðk}, 
F 


这 里 [v^] sz v^ EETA A F RERE, CA F cOn, [vt] = v, Fk: K> KK 为 参考 单 
元 到 一 般 矩 形 单元 的 可 逆 仿 射 变换 。 


3 非 线 性 双 曲 型 积分 微分 方程 及 其 逼近 


考虑 下 面 非 线性 双 曲 型 积分 微分 方程 

ut 一 V.fa(X)VYu + f o(X,t,T,u(X,T))Vu(X,r)dr} + f(u), (X,t) e0 x (6,T], 

u(X,0)-0, wX,0-0, XEN, u(X,t) 20, (X,t) € aN x (0,T], 
(1) 
其 中 XX = (z,y)。 方 程 的 系数 满足 下 述 假定 : 

O 存在 正常 数 ago，a1， 使 得 ao < a(X) <a: 

(Gi) b(X,t,r,u(X,7)), (X, t, r,u(X, )) F f(u) E E u Wi E Lipschitz 连续 条 件 且 均 为 
有 界 函 数 ， 并 具有 本 文 论证 所 需 的 各 阶 有 界 偏 导 数 。 

对 应 于 (1) 的 变 分 形式 为 : Ku: [0,7] 一 HN), EE 


| (un, v) + (a(X)Vu, Vv) + ( fd b(u(r)) Vu(r)dr, Vv) = (F(u), v), Yve H4O), G 
u(0)=0, u(0)=0, 
其 中 
(ug) se ri wvdzdy, b(u(r)) - (X,t,r,u(X,r)), u(0)— w(X,0), (0) u(X,0). 
x (2) 的 非 协 调 有 限 元 半 离 散 格式 为 : cR u^ : [0,7] 一 V*， 使 得 
(3) 


(ub, v) + (a(X)Vu^, Vv) + (fe bur(T) Vu^(r)dr, Vv), = (f(u"),vy, VveV^, 
ur(0) =0, u^(0)- 0， 
HER (Ja = DC dre 容易 证 明 
IM - (35b) 
K 


V^ 上 的 模 。 


第 2 期 在 东洋， 王 慧 敏 : 非 线性 双 曲 型 积分 微分 方程 的 各 向 异性 非 协调 有 限 元 逼近 


4 基本 估计 式 
记 包 = 二 4 一 了 uw， 其 中 五 为 V* 上 的 插值 算 子 。 
引 理 1 图 3 
b(u)e L*(W'»), ue Hl(qQ)nH?(Q), 
d 


b(u) € L*(W?*), ue Hl(0)n H*(Q), 
时 ， 对 任意 的 ve Vih， 则 在 各 向 异性 网 格 下 分 别 有 


à 
[E S Go Gas| < cnfol lame balls x chlrullielo 
K 13K n 


或 
ð 
[E f oo vas < erro( lee lulallolls < chulo. 
K oar On 
3|:gp20) 有 下 面 的 估计 
lwlo + hllwlls < cr?lul2, v we HI(Q) N R? (N), 
(Vw,Vv), =0, (Vw,Vv) =0, VveV*. 
引 理 3 Wa(X), blu) 满足 假设 条 件 ， 对 任意 的 ve V^ 有 
(a(X)Vw, Vv)s + (b(u)Vw, Vv), = O(h?)|ula||v||s. 
证 明 令 
NUM 
a= |K| ] ax), 
则 有 


(a(X)Vw, Vv), — y (aVau， Vo)k 十 D ((a(X) — a)Vw, Vv) y 
K K 


IA 


ch y |wh, vli < ch?lulalvlla. 
K 
同 理 可 证 (b(u)Vw, Vv) € ch?|ulal|vi|n. 3IHAE. 


5 ”收敛 性 分 析 
利用 上 述 引 理 ， 我 们 可 以 得 到 下 面 的 误差 估计 。 
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(4) 


(5) 


(6) 


定理 1 设 w 是 (1) 的 解 且 w, wt € L*(H?)(|] L?(W'99),us E€ D(H?), w^ e V+ 是 半 离 散 


问题 (3) 的 解 ， 则 有 


ju = u^ |l ro q.i) T lu 一 ur? |lzeeqz2) < ch(|ul| Los (x2) 十 lucro (2) 十 [uei r2(82))- 


(8) 
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证 明 记 9 二 届 一 wu， 其 中 09€ Vr。 由 (3) 式 和 (1) 式 可 得 误差 方程 ,对 任意 的 ve€E V* 有 
(ub — ug, v) = —(a(X)Vu^ —a(X)Vu, Vv), «(f b(u(r))Vu(r)dr, vo), 
t (f(u*) — f(u), v) 一 ( E b(u^(r))Vu^ (r)dr, ii 
-he gite Y f. (fru us (9) 
在 (9) hito —60, PAX OSIUEA, di T6(0) = 0(0) = 0， 则 有 
Ilo: + 5 (aC V6, v6), = f Kom n * (a(X) Vw, V&)) dt 
i (f b(u(r))Vw(r)dr, V6.(t)) dt 
= [ ( ri b(v^(r))Vé(r)dr, VO: (t)) àt 
0 
P. d ( T (b(u(r)) — blu” (T)))VInulr)dr, Ve: (t)) dt 
d u^ S u t rd ——Ut s) 
+f (69) - fe),6)at [( (he aede) 
B & Du 人 7) x š 
[ x f (netz PT) ar)6dsdt = Ye (10) 
利用 Halder FR, Young 不 等 式 及 引 理 1 中 (4) 式 估计 Gi(i = 1,2,:… ,8) 
Gi Ge < a f (Ius + |u|2)dt + : f (llel + loh? dt; 
t' t 
G2 +G + G4 = (a(X)Vw, V6), — f (a(X)Vw:, VO) „dt + ( f b(u(r)) Vw(r)dr, VO) 
t t Ü 
2 li ( > be(u(7)) Vw(r)dr, vO(t)) dt — / (b(u)Vw, V6), dt 
-([ we (v^ (7))V(r)dr, V0). + f (b(u*) VO, V8) dt 
h T T. 
$ f ( 7 bí(u^ (r)) VO(r)dr, VOC) dt 


E t 
< ejt + ch? (ue f (Q+ luli)ae) +e f Rat (11) 
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iT lel < c ft (eli + lje«ti2)at, i9llo < lle]. ( 见 文献 [四])， 则 有 


Gs 


( f (b(utr)) — blu” (7) V Isu(r)dr, vb) 一 f ((b(u) — b(u^)) V Inu, VO) , dt 


— fi ( f (bs Cu(r)) — belu” Cr) V Iu(r)àr, vO)) dt 
0 0 h 


I^ 


eje] + ch f ulBdt + c f lea; (12) 


Crte = -Dh X90as+ f Qon a(X) 24 gas) dt 
235 (/ bu(r)) D az) eas 
4 DIM [ (alr) Tar)ol)ds) e f p» b( (u) 30ds)d 


t t 
< eon ee (Ia f aiat) +e f Nolar. (18) 


tr T (a(X)V0, V6), > ao||9||2. £& 结合 上 述 各 估计 式 ， 取 适当 小 的 e， 整 理 得 


t t 
lallà + uel < ch? ( f (Iul + Isl) dt + Ilui?) + c / (lelő + I1) dt 
利用 Gronwall 不等式， 则 有 
t i 
lso + lll. chulo ( f Quel Iuda) ). 
0 


有 引 理 2 及 上 式 ， 并 利用 三 角 不 等 式 ， 定 理 得 证 。 

通过 进一步 分 析 ， 我 们 还 可 以 得 到 下 面 的 超 逼 近 和 整体 超收 和 敛 结果 。 

定理 2 设 久 是 (1) 的 解 且 w e L?(H?), u, e D(H*)(]L*(W?9?), un e DP(H?), w^ € 
V* 足 半 离 散 问题 (3) 的 解 ， 则 有 


lu” — Tawra) S eh? (ullos) + lules) + furl), (4) 
|ru” 一 u|| r-» (U-I4) < ch? (||u]| res r2) T lus r2 3) 十 |utzlzaz(z2))， (15) 


IEP nu 29 u WARIH., Do, 为 利用 文献 [6] 中 的 思想 构造 的 插值 后 处 理 算 子 。 

证 明 采用 定理 1 的 证 明 方 法 ， c SA Gs 利用 引 理 3，Gz, Gg 利用 引 理 1 中 
的 (5) 式 ， 则 (14) 式 得 证 。 利 用 文献 [6] 中 的 思想 ， 类 似 于 文献 构造 插值 后 处 理 算 子 Dos, 
则 (15) 式 得 证 。 从 而 定理 得 证 。 
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The Nonconforming Finite Element Approximation to the Nonlinear 
Hyperbolic Integro-differential Equation on Anisotropic Meshes 
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Abstract: The approximation of a rectangular nonconforming finite element method for a kind of 
nonlinear hyperbolic integro-differential equations is discussed under the anisotropic meshes, the con- 
vergence analysis and the error estimate of finite element solution are presented for the semi-discrete 
Scheme. The superclose property is derived through the new technique and sharp estimates when the 
exact solution is appropriately smooth. At the same time, based on the interpolated postprocessing 
trick, the global superconvergence is obtained. The results of this paper show that the regular condition 
on the meshes and the dependence on the Ritz-Volterra projection in traditional finite element analysis 
are not necessary, and thus extend the applications of nonconforming finite element methods. 
Keywords: hyperbolic integro-differential equation; anisotropic meshes; nonconforming finite element; 
semi-discrete; superconvergence 
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